Pure spin current in graphene NS structures 
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We demonstrate theoretically the possibility of producing a pure spin current in graphene by 
filtering the charge from a spin-polarized electric current. To achieve this effect, which is based on 
the recently predicted property of specular Andreev reflection in graphene, we propose two possible 
device structures containing normal-superconductor (NS) junctions. 
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I. INTRODUCTION 

The recent experimental realization of conducting two 
dimensional monolayers of graphite [H, 0, [H 3, a l so 
known as graphene, offers the promise for new electronic 
devices. One conceivable use for graphene is in spintron- 
ics @, where the lack of nuclear spin interaction ( 12 C 
has no nuclear spin) could offer the ability to maintain 
spin coherence over larger distances than in conventional 
semiconductors. For progress towards this goal, it is es- 
sential to have simple and reliable means to transport 
spin in graphene. 

We address this issue by proposing a prescription for 
producing a pure spin current in ballistic bulk graphene. 
Spin currents have already been predicted to arise in 
graphene due to spin-orbit coupling Q and the Quan- 
tum Hall Effect [7]. In both cases, the spin currents are 
due to counter-propagating edge states of opposite spin. 
Our proposal makes use of the recently predicted specu- 
larity [|[ of Andreev reflection [§| in graphene to produce 
a pure spin current in bulk. This is accomplished us- 
ing structures containing normal-superconducting (NS) 
boundaries, which filter the charge out of a current 
of spin-polarized quasiparticles, leaving behind a pure 
spin current. To be specific, we consider two device 
paradigms: 1) a V-junction geometry with opening angle 
appropriately tuned, 2) a channel with a superconductor 
at one boundary and a normal edge at the other. An ad- 
vantage of such devices is the large number of transmit- 
ting channels in the bulk, offering the possibility of rapid 
spin accumulation. Also, our proposal does not require a 
magnetic field, nor does it rely on the spin-orbit gap. In 
the present case, however, it is crucial to first generate 
a spin-polarized electric current, which could conceivably 
be done by contacting the system to a ferromagnetic lead, 
as has been done with carbon nanotubes [lfj. 

Our description of proximity effects in graphene fol- 
lows that of Ref. Q. Later publications, Refs. [Til H3|. 
used this approach to discuss the Andreev spectrum and 
Josephson effect in NS (SNS) structures in graphene. A 
very recent paper [13| makes explicit use of the specu- 
larity of Andreev reflection in graphene by considering 
the neutral excitations propagating along a narrow SNS 
channel. These authors propose a device that is similar to 




FIG. 1: V-junction with NS interface. Polarized electrons 
are injected through the lower arm and the same polarization 
is transferred to the upper arm by the electron-hole beam. 

the one considered below in Sec. IIV[ but suggest investi- 
gating the thermoelectric effect to observe the chargeless 
excitations. In contrast, we analyze the spin transport. 
Our approach also lends itself to an analysis of the de- 
viations from perfect charge filtering in a channel, which 
is done in Sec. [V] Finally, an elegant method to pro- 
duce pure spin currents in conventional semiconductors 
was suggested in Ref. where the spatial separation 
of electron and hole trajectories is caused by tunneling 
through a superconductor. 



II. ANDREEV REFLECTION IN GRAPHENE 

The electron wave function in graphene is described 
by a two component (pseudo-)spinor ip. Its spin-up and 
spin-down components correspond to the quantum me- 
chanical amplitudes of finding the particle on one of the 
two sublattices of the honeycomb lattice. The low energy 
physics of graphene is governed by two so-called Dirac 
points in the spectrum, located at the two inequivalent 
corners K, K' of the Brillouin zone. The spinor wave 
function for low energy excitations in (lightly-doped) 
graphene decomposes into a sum of two waves oscillat- 
ing with different wave vectors tp — e lKr (t)+ + e lK ' r a L >_. 
The smooth envelope functions cf>± satisfy the two- 
dimensional Dirac equation (l5j described by the Hamil- 
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tonian 



described by a 2 x 2 reflection matrix, R A , 



c(a x p x ± cr^) + U(x), 



(1) 



where c is the Fermi velocity and p = — iftV. In regions of 
constant U this equation defines a conical energy band, 
or valley, e — U = ±c|p|. The Pauli matrices o~ XjV per- 
mute electrons between two triangular sublattices of the 
honeycomb lattice. The two signs in Eq. fTJ) (+) and (— ) 
correspond to the two valleys K and K' . 

We consider a graphene sheet in the x — y plane, with 
the region x < covered by a conventional supercon- 
ductor (Fig. [T]). Following Ref. § we assume that in 
this configuration a pair potential A(x) — AqQ(-x) can 
be induced [l6| in the graphene sheet by the proximity 
effect, accompanied by a sufficiently strong shift in the 
scalar potential U(x) — —UoQ(—x), so that Uq ^> A . 
The reflection at the NS interface (x = 0) is described by 
a separate four- dimensional Dirac-Bogoliubov-de Genncs 
equation [8( for each valley 



H ± - E F A(x) 
A(x) E F - H ± 



(2) 



where the two spinors u and v represent the electron 
and hole components of the <p± wavefunctions for H±, 
respectively. 

In addition to carrying pseudo-spin, electron-hole, and 
valley indices, the quasiparticle wave function should also 
describe the usual spin. Since for e < A no spin can be 
injected into the superconductor, the spin of the incident 
electron is transferred to the reflected particle. 

For several decades it was considered a basic feature 
of Andreev reflection [9( that the hole produced upon 
electron-hole conversion retraces the incident electron's 
trajectory. Even in the traditional materials, however, 
this repetition of trajectories is exact only for zero exci- 
tation energy e in Eq. @ . At finite excitation energy the 
two trajectories do not quite coincide, leading to inter- 
esting effects in semiclassical dynamics [l?} and the spec- 
trum [l8| of Andreev billiards. In graphene, the Andreev 
reflection described by Eq. @ has the standard form 
only if e <C Ep . Because of the possibility to tune the 
Fermi energy close to the Dirac point, one may reach in 
graphene the regime e > so that an incident electron 
in the conduction band produces an Andreev reflected 
hole in the valence band (opposite side of the cone in the 
Dirac excitation spectrum). This valence band hole has 
the same velocity along the NS interface as the incident 
electron, and consequently is specularly reflected [8| . For 
vanishing Fermi energy, both the electron and hole com- 
ponents of the reflected wave follow precisely the same 
trajectory. We now restrict ourselves to this most inter- 
esting case of Ep = 0. 

The transformation of a quantum superposition of in- 
cident electron and hole into the corresponding super- 
position of outgoing waves upon Andreev reflection is 
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With the notation a — arctan \p y /p x \, e = c^Jpl, + p 2 yl 
and £ = \J A^ — e 2 one has [HI 
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£ + i£ cos a 



(4) 



III. 



NS INTERFACE IN V-JUNCTION 
GEOMETRY 



Characteristic of specular Andreev reflection in 
graphene is the spatial separation of incident and re- 
flected electron-hole beams. The simplest device that 
makes use of this property is a V-junction, as shown in 
Fig. 1. Suppose one can inject a spin-polarized collimated 
monoenergetic beam of electrons through one arm of the 
junction. [Monoenergetic beams may potentially be pro- 
duced by resonant transport through a graphene quan- 
tum dot (QD) as was done for semiconductor QD's 
e.g. in Ref. [2l|.] The injected electron will be either nor- 
mally or Andreev reflected at the NS interface and the 
created quasiparticle will escape through the second arm. 
Since no spin can penetrate through the superconductor, 
all the polarization of the injected beam is transferred to 
the second arm. On the other hand, the average charge 
of the quasiparticles reflected towards the second arm 
depends on the excitation energy. One obtains zero re- 
flected total charge by setting \r\ = \r A \- This constrains 
the incident particle energy to be 



e = A cot a. 



(5) 



Alternatively, one may consider the injection of a col- 
limated beam of electrons with all possible energies in 
the range < e < eV. For Dirac particles and a fixed 
angular spread of the beam, cancellation of the reflected 
charge requires 
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\r A \ ede. 
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One may think here about electrons leaving the bi- 
ased (eV) graphene half-plane through a narrow slit and 
then both angle- and width-collimated by a second slit. 
Substitution of Eq. (HJ) into Eq. ^ leads to the transcen- 
dental equation 
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with the solution 



eV = 1.59 An cot a. 



(7) 



(8) 
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Experimental realization of pure spin current in a 
graphene V-junction is limited by the requirement of a 
collimated electron beam. This difficulty is avoided auto- 
matically in the setup proposed below, where the filtering 
of charge current takes place because of multiple Andreev 
reflections in a NS channel. 



IV. CHARGE FILTERING IN A LONG NS 
CHANNEL 

We consider a normal graphene strip of width W and 
length L 3> W, which is formed by covering the region 
x < of a graphene half-plane — oo < x < W with 
a superconductor (see Fig. 2). Spin-polarized electrons 
are injected at y — 0. We restrict our discussion to the 
geometric optics limit, A <C W, where the transmission 
through the channel may be considered in terms of in- 
dividual particle trajectories. Here A = 2-KficjeV is the 
de Broglie wavelength for Dirac electrons with energy eV. 
In addition we consider the low- voltage limit eV <C Ao, 
since in this limit the normal reflection at the NS inter- 
face is suppressed (see Eq. Q) 



\r/r A \ 2 ~ (eV/Ao) 2 « 1, 



(9) 
(hole- 



and we may assume a perfect electron-hole 
electron) conversion by Andreev reflection. 

The number of transmitting channels N(e) for the en- 
ergy e is determined by the width of the strip W and the 
range of variation of transverse momentum \p x \ < e/c, 



eW 



(10) 



Here the first factor of 2 accounts for the two valleys in 
graphene. We do not add another 2 for spin, since we 
assume injection of polarized current. [In the limit given 
in Eq. ([9]), the possible inter- valley mixing upon charge- 
and spin-preserving reflection at the normal edge, see e.g. 
Refs. [iH [23j], is not important for our calculation.] As 
usual, each open channel adds Go = e 2 /h to the differ- 
ential conductance. 

The total current injected into the strip is [2~j ] 



N(eV')G dV' 



< 2 tV-]V _ ^ e 3 V 2 W 



h chir 



ch 2 



Associated with each incoming electron is a certain quasi- 
particle trajectory in the channel as depicted in Fig. 2. 
In the low- voltage limit, Eq. ©, the outgoing particle is 
an electron if the number of reflections from the NS inter- 
face is even, and a hole if the number of reflections is odd. 
The crucial point is that for a long channel the contribu- 
tions to the total current of these two kinds of trajectories 
effectively equilibrate, so that the total charge transfer 
vanishes after averaging over initial angles. While the 
charge current upon such equilibration escapes to the su- 
perconductor, the spin current continues to flow along 





FIG. 2: NS-strip. Charge current injected from contact 1 
escapes to the superconductor 3, leaving behind a pure spin 
current flowing towards electrode 2. 
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FIG. 3: The average charge per quasiparticle after passing 
through the strip (length L, width W , L/W = 10) as a func- 
tion of transverse momentum p x . Very close to \p x \ = e/c 
the function Q(jp x ) could not be shown because of the fast 
oscillations. 



the strip. For a 100% polarized injected beam the rate 
of spin transfer is simply found as 



dS_ 

dt 



hi 
2 e' 



(12) 



V. SOURCES OF INCOMPLETE CHARGE 
FILTERING 

For a channel of finite length the numbers (weights) of 
trajectories with even and odd numbers of Andreev re- 
flections do not quite coincide, leading to a finite charge 
current. In this section we first determine the fraction 
of charge per quasiparticle (Q) remaining in the beam 
after passage through the NS channel, due to this purely 
geometrical effect. Since the normal boundary acts as a 
perfect mirror, we may effectively consider a strip of dou- 
bled width, < x < 2W, having NS-interfaces on both 
sides. The fluxes of electrons injected into this doubled 
strip with p x < and p x > are physically equiva- 
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lent. The charge per quasiparticle with incident angle a 
changes linearly along the channel from — e at y = to 
+e at 



y = yo = 2VFtana. 



(13) 



At y = y all electrons injected at angle a are converted 
to holes. At yo < y < 2yo the holes are converted back 
to electrons. Charge here changes back (linearly) from 
+e to — e, and so on. 

Figure 3 shows the average charge per quasiparticle 
for electrons injected at a fixed angle a at various initial 
positions, < x < 2W, as a function of p x c/e — cos a for 
L = 10W. The segments of the function Q(p x ) are given 
by (n = 0,1,2,...) 



Q{p x ) = e(-l)« 



L 



C\Px\ 



w ^ 2 - (c Px y 



- (271 + 1) 



(14) 



Averaging (integrating) further over p x and introducing 
a new variable r = cp x / \/ e 2 — (cp x ) 2 , we find the charge 
fraction (Q) per quasiparticle remaining in the beam to 
be 



(Q) 



F | — T 

W 



dr 



(1 + r 2 ) 3 / 2 



(15) 



Here we have introduced the sawtooth function F(x + 
4) = F(x) and F(x) = \x\ — 1 for \x\ < 2. Keeping only 
the first harmonic in the Fourier transform of F(x) in the 
integral in Eq. (fT5|) gives 



(g) "-VV| eXP V 2TIV 



Ltt\ 



(16) 



Remarkably, the corrections to charge neutrality are ex- 
ponentially small in the limit L 3> W. Since this in- 
complete cancellation of charge is of geometric origin, 
Eq. (fl6|) does not depend on the voltage bias eV. 

Another source of incomplete charge filtering in the NS 
channel is the coexistence of Andreev and normal reflec- 
tion for finite quasiparticle energies (e ~ eV < Aq). To 
find the charge transfer in this case we add to the reflec- 
tion matrix, Eq. ([3]), a part describing particle propaga- 
tion in graphene between Andreev reflections 



R — RqRa- 



(17) 



The form of Rq is sensitive to the details of quasiparticle 
reflection at the normal side of the strip. Depending on 
the microscopic structure of the graphene edge, reflection 
from it may or may not introduce transitions between the 
two valleys K, K'. Here we consider only the latter case 
as an example. This also means that Rq is a 2 x 2 diagonal 
matrix. 

In the case of decoupled valleys the boundary condi- 
tions may only have the form (see Ref. [25] for the general 
situation) 



(aipi + bijj 2 )\ x =w = , a, b = const 



(18) 



where ip\ and ip2 are the up and down components of ei- 
ther the particle (u), or hole (v) wave function. Particle 
number conservation imposes certain restrictions on the 
allowed values of a and b. Below we consider two exam- 
ples of such boundary conditions, demonstrating both 
existence and absence of corrections to charge filtering 
due to finite voltage. 

1. The first option is a/b — ±i (for the edge along 
the y axis). This boundary is realized if the particle 
confinement is achieved by adding a term with large mass 
a z Mc 2 to the Dirac equation ([!]) at x > W [29. |27|. 
Straightforward calculation for such a boundary gives 



Ro = ze^^diagtl, -1) , (R Ra) 2 oc I. 



(19) 



The particle-hole superposition produced upon Andreev 
reflection of an electron returns to the pure electron 
state after the second reflection from the superconduc- 
tor. The expectation value of the quasiparticle charge 
along a given trajectory according to Eq. (Till|) switches 
from — e to (\r A \ 2 — \r\ 2 )e and back after each reflection 
from the NS-interface. [The average charge of the beam 
with initial angle a now changes linearly between these 
two values, not between — e and e as in Eqs. (| 141 1 5[) .] Av- 
eraging over angles and energy using the exact Eq. (j4|) 
gives the charge per quasiparticle transmitted through 
the graphene channel for eV <§C Aq as 



(Q) = 



7r eV 
3 AT 



(20) 



We note that according to Eq. (fj| the condition Ir/r^j <C 
1 is violated for grazing trajectories, having tt — a ~ 
s/Aq (j4|). The small statistical weight of these grazing 
trajectories is the source of the small value of (Q). 

The result Eq. ([20]) was found in the limit L > W. 
However, further increase of the channel length does not 
lead here to charge relaxation, in contrast with the ge- 
ometric correction, Eq. (|16|) . The possibility for finite 
charge current, Eq. (f20|) . to flow without leaking along a 
(arbitrarily) long NS-interface is very counterintuitive. 

2. The second possibility consistent with particle con- 
servation Eq. (fTS]) is a x b = (i.e. either a — 0, or b = 0). 
Boundary conditions of this type describe the bulk enve- 
lope functions in graphene with a zigzag edge Q. Re- 
flection from the normal edge is now described by a pure 
phase matrix Rq oc diag(l, 1). Therefore the product of 
n reflections reduces to 



(R) n (x {RaT oc I cos n(3 + ia x sin n/3, 



(21) 



where cot/3 = — etana/Ao and the matrix o~ x inter- 
changes particles and holes. Eq. (f2Tj) leads to a uni- 
form (to exponential accuracy, as in Eq. (|16[1 ) mixing 
of particles and holes after many reflections at the chan- 
nel boundaries for any value of (3, i.e. at any e < Aq. In 
contrast to the previous example, the charge per quasi- 
particle transmitted through the graphene channel for 
L >> W is zero even at finite values of e/ Aq. 
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VI. CONCLUSIONS 



In summary, we have proposed two basic device 
concepts for conversion of polarized electric current 
into pure spin current in ballistic bulk graphene. The 
second device, where the filtering of charge originates 
simply from the equilibration of the number of trajec- 
tories experiencing an even or odd number of Andreev 
reflections, seems especially promising. We expect this 
idea to be easily generalized beyond ballistic transport 
and for different geometries. We stress that neither 
proposal requires a 100% polarization of the incident 
beam, although the spin transfer rate will of course 
depend on the initial polarization. 
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